Shear-driven solidification of dilute colloidal suspensions 
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We show that the shear-induced solidification of dilute charge-stabilized (DLVO) colloids is due 
to the interplay between the shear-induced formation and breakage of large non-Brownian clusters. 
While their size is limited by breakage, their number density increases with the shearing-time. 
Upon flow cessation, the dense packing of clusters interconnects into a rigid state by means of 
grainy bonds, each involving a large number of primary colloidal bonds. The emerging picture of 
shear-driven solidification in dilute colloidal suspensions combines the gelation of Brownian systems 
with the jamming of athermal systems. 
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Shear-driven solidification of diluted colloidal suspen- 
sions has dramatic impact on their applications, from 
industrial polymer production to natural microfluidic de- 
vices 0, IH , and is a prototype of non-equilibrium transi- 
tions. If intcrparticle interactions are purely attractive, 
the applied shear stress may break down aggregates and 
fluidize the material 0, H|. However, the colloidal par- 
ticles are often stabilized by electrostatics [l|, with no 
tendency to aggregate at rest, and high shear rates may 
ultimately promote aggregation in competition with the 
electrostatics. The interplay between these two tenden- 
cies may lead to persistent structures [H and this shear- 
induced aggregation might be as dramatic as a complete 
solidification of even diluted suspensions, hence seriously 
affecting the material and rheological properties. Al- 
though this is a widely reported phenomenon in both 
artificial and living systems [H, Q, there is little un- 
derstanding of the solidification mechanism. This is due 
to the difficulty to monitor the system with real-space 
optics or scattering techniques at high shear rates d, Q. 
To overcome these obstacles, we have designed an ex- 
perimental protocol exploiting charge-stabilized colloidal 
particles which interact via a typical DLVO potential [l| 
where an energy barrier coexists with a deep, shorter- 
ranged attractive well. In our setup, shear-induced ag- 
gregation can be monitored in a fully controlled way, al- 
lowing for the first time to rationalize the effect of the 
shear stress from an initially dilute suspension to the fi- 
nal solid. By combining light scattering, rheology, and 
microscopy data we formulate a model for the solidifi- 
cation mechanism. The effective packing fraction of the 
aggregates formed under shear increases with time. Upon 
flow cessation their dense packing is progressively frozen 
into a rigid structure by the formation of grainy (i.e., 
multiple) colloidal bonds which are responsible for the 
fairly high shear moduli observed. 

Experiments. The system consists of colloidal particles 
at a fixed colloid volume fraction of cf> ~ 0.21. For the 
effective DLVO interaction [l|, the attractive well depth 
is ~ 40ksT and the repulsive barrier ~ 60/sbT. The col- 



loidal particles are surfactant-free polystyrene-acrylate 
latex spheres charge-stabilized by the charged groups of 
the initiator 0] ■ The nearly monodisperse particles have 
mean radius a ~ 60nm as determined from both dynamic 
and static light scattering (using a BI-200 SM goniometer 
system, Brookhaven Instruments, NY). YlmM of NaCl 
were added to weakly screen the electrostatic interaction 
barrier, characterized by Zeta-potential measurements. 
A strain-controlled ARES rheometer (Advanced Rhco- 
metric Expansion System, TA Instruments, Germany) 
with Couette geometry has been employed to induce the 
shear flow under shear-rate control and to measure the 
rheological properties. We have initially sheared the sys- 
tem for a varying time t\ at a fixed shear-rate (1700s -1 ). 
For each n, upon flow cessation, we sampled the shear 
cell and analyzed the system by laser light scattering 
(LLS) using a small-angle light scattering Mastersizer 
2000 instrument (Malvern, U.K). The LLS analysis can 
be done off-line since the aggregation docs not evolve 
on the time scale of the analysis in the absence of shear. 
The samples were diluted at 7 = to such an extent that 
multiple light scattering does not affect the measurement. 
After each t\ we have also performed rheological shear- 
sweep and frequency-sweep tests. Since the height of the 
energy barrier for bond-breaking is ~ lOO/c^T, the clus- 
ters are mechanically stable even under the fairly high 
shear rate of the Couette cell, hence it is clear that they 
cannot de-aggregate during the off-line analysis. From 
the scattered intensity I(q) we obtain the average struc- 
ture factor of the aggregates (S(q)} 3 ( n °t to be confused 
with the one of the whole suspension) present in the sys- 
tem at the time of sampling. 

Results and discussion. The starting point of our anal- 
ysis is the time evolution of the systems sheared at a 
constant shear rate 7 = 1700s . The viscosity shows 
a sharp upturn in time, after an induction delay due to 
the activation barrier in the microscopic aggregation ki- 
netics between two colloidal particles [1| . The barrier for 
shear- induced aggregation, in fact, decreases exponen- 
tially upon increasing the shear rate and vanishes upon 
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FIG. 1. Main frame: Time evolution of the viscosity r\ at 
7 = 1700s -1 an tf> ~ 0.21. Insets: left - rj/r]o as a function 
of cj> e f f from experiments (circles) and the extended Einstein 
formula (line); right - Optical micrograph upon flow cessation. 



reaching a certain cluster size above which the aggre- 
gation kinetics is much faster and goes with the cube 
of the cluster size [lj]. Hence a kind of self-accelerated 
kinetics, leading to larger and larger aggregates, is ex- 
pected once clusters of a critical size are formed. On 
the other hand, the action of shear is also well known to 
cause fragmentation, which should ultimately lead to a 
maximal cluster size, dictated by the mechanical balance 
between the stress imposed by the flow on the cluster and 
its mechanical response. Nevertheless we observe (FigfTJ 
main frame) that the viscosity continues to increase until 
the instrument stress overload threshold is reached. We 
have stopped the shearing after different time duration 
Ti (pre-shearing times) . For each n we have used optical 
microscopy and LLS to investigate the structure of the 
resulting suspension. An optical micrograph upon flow 
cessation at the largest t\ is shown in the inset (right) 
of Fig|T] From the small angle light scattering analy- 
sis [10j, it is evident that already after a short time the 
size distribution of aggregates is strongly bimodal (i.e., 
primary particles coexist with large aggregates). The 
Guinier plot gives a typical aggregate radius, resulting 
from the competition between aggregation and breakage, 
R g ~ 35±3/zm which remains constant with time [loj ]. in- 
dicating that the aggregation under shear rapidly leads to 
an optimal cluster size. From the power-law regime of the 
scattering curves we extract the fractal dimension df of 
the large clusters: as a function of n [Io| df also rapidly 
reaches a plateau at ~ 2.7±0.1, a fairly high value which 
is typical of shear-driven aggregation [6( where breakage 
events and restructuring induced by flow stresses com- 
bine JjJ. We have measured the fraction x(t) of primary 
colloidal particles converted to clusters and found that it 
steeply increases [loj ]. The emerging picture is that un- 
der shear the system is constantly generating new clusters 




FIG. 2. (color online) Shear-sweep curves obtained after pre- 
shearing at 7 = 1700s -1 for n = 4054s, 7420s, 7500s, 7520s 
and 7555s (from bottom to top). 



with approximately the same size R g and fractal dimcsion 
df. The fast increase in the number density of clusters 
which arc fractal and porous is associated with a rapid 
increase of their effective packing fraction: this could be 
the reason for the sharp non-linear increase of the vis- 
cosity in Fig|TJ To test this hypothesis, we have esti- 
mated the clusters effective packing fraction 4> e f f through 
the relation <f> e ff(t) = x{t)4'k~ 1 {Rg/a) 3 ~ df , where k is 
a geometric prefactor close to one and (j) is the initial 
volume fraction of primary particles In the inset 
(left) of FigJT] we plot the viscosity rj(t) (circles), nor- 
malized by t]q, as a function of 4> e ff(t). We have also 
calculated the high shear viscosity of an equivalent sus- 
pension of hard spheres of the same linear size at the 
same volume fraction cf> e ff(t), using the Einstein for- 
mula properly extended to high 4> e ff(t) [HI which gives 
77 ~ ?7o exp (5(/) e //(i)/2) (the full line in the same inset). 
The agreement indicates that the increase of the viscos- 
ity under shear in FiglT] can be ascribed indeed to the in- 
creasing packing fraction of the clusters, which hydrody- 
namically behave as hard spheres due to the fairly high d / 
flit ) : the initial dilute colloidal suspension has changed, 
under shear, into a suspension of non-Brownian aggre- 
gates whose packing fraction increases with the shearing 
time. 

After each of the pre-shearing times, we also perform 
a shear-sweep experiment where 7 is varied but kept be- 
low 7 = 1700s" 1 . This guarantees that the aggregates 
formed during the pre-shearing are mechanically stable 
and do not break up during the shear-sweep [13, 
The data are plotted in Figj2] The curves correspond 
to ti = 4054s, 7420s, 7500s, 7520 and 7520s from bot- 
tom to top, and to 4>eff(i~i) increasing accordingly. For 
the shortest n, a modest shear-thinning is followed by 
a Newtonian plateau. Upon increasing n, the viscosity 
curves are shifted to higher values and a steep shear- 
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FIG. 3. (color online) The measured stress as a function of 
shear rate (symbols) and the extended Hershel-Bulkley curves 
(lines) 
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FIG. 4. (color online) G" (filled symbols) and G" (open 
symbols) measured after ri = 4054s, 6350s, 7400s and 7480s 
(from bottom to top) as a function of the frequency ui in units 
of the diffusion time in the dilute limit, to = a 2 /Dq. 



thickening appears: the rheological behavior of our sys- 
tem upon increasing n is qualitatively the same as the 
one of dense non-Brownian suspensions upon increasing 
volume fraction [l9|, [2(| • Note that flow inhomogeneity 
may occur in this type of measurements 2JJ, l22| and we 
cannot exclude its presence here (although MRI studies 
in a very similar system [23| detected shear banding only 
for 7 < 60s -1 ). Thus Fig(2]gives only macroscopic rheo- 
logical information, possibly dependent on the geometry. 
A shear thinning followed by a shear thickening regime 
has also been observed in some gelling colloidal suspen- 
sions [24[, where it has been explained in terms of an 
increase of the clusters number density due to break-up. 
Here we have made sure that the shear-sweep is done 
at 7 lower than the value needed for breaking up the 
clusters, therefore we rather associate the shear thick- 
ening to the formation of inter-clusters bonds, possibly 
also promoted by hydrodynamic interactions [25j and/or 
from a residual inter-cluster attraction: the shear rate 
is too low to break the pre-formed aggregates but can 
promote, instead, the formation of further bonds among 
them. At the largest <t> e ff (i- e -> r i) the shear-thickening 
is no longer visible. The behavior of the stress intensity 
a as a function of 7 plotted in FigJ3] sheds some light 
on this rheological behavior. For large 4> e ff the curves 
can be well described with an extended Hcrschcl-Bulklcy 
typ e of behavior, which accounts for shear-thickening 
[261 ] . a(j) = o y + 017 1 / 2 + a27 1 / e where a y is the yield 
stress, ai and 02 are parameters depending on 4> e ff and 
e < 0.1 is typically obtained for the steep shear thickening 
of densely packed non-Brownian suspensions [26| . This 
indicates that the onset of the yield stress occurs in the 
range of effective volume fractions 0.50 < cf) e ff < 0.56. 
At (f) e ff ~ 0.66, i.e. the largest value obtained here, the 
yield stress has now a value close to the upper stress 



limit of the shear thickening and this makes the shear 
thickening no longer visible. The same phenomenon has 
been observed indeed in densely packed non-Brownian 
suspensions. Hence, in a dilute colloidal system at a solid 
fraction <fi = 0.21 we get the whole rich phenomenology 
observed in non-Brownian suspensions upon varying the 
solid fraction in a broad range. 

The emergence of a yield stress at sufficiently high <p e f / 
raises the question of how a fully solid state may form. 
Upon flow cessation after each pre-shearing time t\, we 
have also performed dynamic frequency sweeps. The re- 
sults are reported in FigJU At the lowest T\ considered 
(which is however close to the point where the high-shear 
viscosity starts to rise dramatically) the system still be- 
haves liquid-like as the loss modulus G" is always slightly 
larger than the elastic modulus G' . Upon increasing t\ 
and hence <f> e ff, G' takes over with respect to G" and 
the ratio G'/G" becomes larger with n, while remaining 
mostly constant with the frequency. This is a behav- 
ior typically observed in colloidal gels [H, 28, [29[ and 



suggests that, for the longest T\, upon flow cessation the 
clusters are not only densely packed but can also connect 
into a spanning, stress-bearing structure (3pj . During the 
shearing, in fact, the aggregates are maintained in a fluid 
state by the imposed high shear rate, whereas upon ces- 
sation of flow the system is subjected to a rapidly de- 
creasing stress until a quiescent state is reached. Hence 
inter-aggregates connections can gradually form and be- 
come permanent as the zero stress is reached: provided 
that <p e f f is high enough, a cohesive solid random packing 
of aggregates will form. To test this picture, we have cal- 
culated the elastic modulus of this disordered solid using 
the approach derived in (3f^ . giving a quantitative esti- 
mate of the contribution to elasticity of inter-aggregates 
connections for an amorphous close-packed assembly of 
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aggregates of cohesive particles: 



G' = (2/5n)K</) eff zR 



(1) 



where z is the aggregate average coordination number, 
R g is their linear size and k is the stiffness of the inter- 
aggregate connections formed upon flow cessation. We 
use this result to estimate the elastic modulus of our ag- 
gregate gel formed upon flow cessation at (f) e ff = 0.66, 
corresponding to the largest effective cluster packing frac- 
tion and the most rigid state observed, where the affinc 
assumption underlying Eq. dTJ is reasonably applicable. 
The irregular aggregates morphology, visible also by op- 
tical microscopy (inset of FigQ}, indicates that each con- 
nection is actually composed by many colloidal bonds 
(grainy contacts). In order to estimate the number of 
such bonds, we consider that for a fractal aggregate of 
radius R g and fractal dimension df, the total number of 
particles is N c = (R g /a) df , a being the particle radius. 
Hence a(dN c /dR g ) = df/aNc- df 1 ^ df gives the number 
of particles added to the outermost layer of an aggre- 
gate of radius R g . Using a ~ 60tito, R g ~ 35 ± 3/xm 
and df ~ 2.7 for our system, we get the number of 
particles on the surface of the aggregates ~ 136 • 10 3 . 
With an average number of 7 grainy contacts per aggre- 
gate, as for densely packed spheroidal objects |3lf . we 
can estimate that each of them involves n ~ 19.4 ■ 10 3 
particles, consistent with the value n ~ 21 • 10 3 ob- 
tained measuring the contact area from optical micro- 
graphs. On this basis, we estimate k = nn, where 
K = (d 2 U DLVO /dr 2 ) r=rmin ~ 2 • 10- 5 N/m. Hence we 
finally obtain for the elastic modulus of our dense gel 
formed upon flow cessation, G ~ 760Pa, consistent with 
the value measured experimentally G" ~ 843Pa at the 
largest <j) e ff (see Fig[4]). 

Conclusions. Our experiments rationalize the shear- 
induced solidification of a dilute, stable colloidal sus- 
pension for the first time in a fully controlled way. 
Large aggregates of a typical size are continuously gen- 
erated under shear and behave hydrodynamically like 
non-Brownian hard-spheroids. Varying the shearing 
time leads to the same rheological response of dense 
non-Brownian suspensions upon varying the solid frac- 
tion. Upon flow cessation, these aggregates can even- 
tually form cohesive random packings where each inter- 
aggregate bond involve a large number of colloidal bonds. 
Such solidification mechanism is thus a hybrid between 
colloidal gelation [32| and the packing-driven jamming 
[33| of non-Brownian suspensions (pastes, slurries). This 
scenario gives a novel insight into the complexity of 
the shear-induced solidification of colloidal dispersions of 
practical relevance. 
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LASER LIGHT SCATTERING (LLS) ANALYSIS 

After each pre-shearing, made at constant shear rate 
7 = 1700s -1 for a duration r%, we used LLS to investigate 
the structure of the suspension in terms of the structure 
factor which is shown in FigJS] Immediately after sam- 
pling the rheometer, the samples were diluted to such 
an extent that multiple light scattering does not affect 
the measurement and such that the aggregates behave 
as nonintcracting scattcrers. At the shortest t\ consid- 
ered (i.e., T\ < 1000s), the aggregate structure factor is 
flat, i.e. the amount of aggregates is negligible 0. Upon 
increasing n, the structure factor develops a plateau at 
low wave- vector q, indicative of large clusters, followed 
by the power-law regime of fractal objects |3[. Up to 
T\ = 4893s, the power-law regime is followed by a fiat 
tail which indicates that the system is still dominated 
by free primary colloidal spheres coexisting with a few 
very large clusters and negligible amounts of clusters of 
intermediate size. Confirmation of this bimodal size dis- 
tribution of the aggregates comes from the analysis of the 
samples after filtration with a 5/xm-pore size filter which 
allows us to remove the large clusters from the sample. 
The filtered sample is then analyzed by LLS and found to 
be composed only by oligomers (monomers, dimers, and 
trimers) 0, H| . The gyration radius of the large clusters 
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FIG. 6. (color online). Main frame: Average radius of gyra- 
tion R g of the clusters (from the Guinier plot) as a function 
of n. Inset (left): Average fractal dimension of the clusters 
df as a function of t\. 



still applicable, as the lowest scattering angle of our in- 
strument is 0.0145°, i.e. well below the limit of 0.1° corre- 
sponding to the typical value of R g . From the power-law 
regime of the scattering curves we also evaluate the frac- 
tal dimension df of the large clusters as a function of 
the shearing time (inset of FiglBJ). The time evolution 
of df also rapidly reaches a plateau at df ~ 2.7. Such a 
fairly high fractal dimension is typical of the shear-driven 
aggregation mechanism [f| and also results from the re- 
structuring induced by flow stresses and breakage events 
Q. Finally we have measured the fraction \{t) of pri- 




FIG. 5. (color online). The scattered intensity of the dilute 
system, I(q), divided by the independently measured form 
factor of the latex spheres, P(q), yields the average structure 
factor of the aggregates < S(q) >= I(q)/P(q) [f|, for increas- 
ing pre-shearing time n (from bottom to top). 

has been estimated from the scattering curves at different 
shearing times through the Guinier plot and is reported 
in the main frame of Fig|6] The observed plateau in 
the cluster linear size clearly indicates that the aggrega- 
tion under shear rapidly leads to a maximal cluster size 

R g ~ 35 ± 3fim. The Guinier analysis in this regime is mary colloidal particles converted to clusters and found 
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FIG. 7. (color online). Fraction \ of primary colloidal par- 
ticles converted to large aggregates at constant shear rate 
(7 = 1700s -1 ) as a function of the pre-shearing time n. 
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that it steeply increases with the pre-seharing time as 
shown in Fig|7] The fact that the size and morphol- 
ogy of the clusters remain constant with time indicates 
that under shear the system is constantly generating new 
clusters with approximately the same size (~ 35 ± 5fim 
of radius) and df ~ 2.7±0.1. The fast increase in the 
number density of clusters which are fractal and porous is 
associated with a rapid increase of their effective packing 
fraction. 
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